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ABSTRACT 

An extension of the Tridiagonal Reduction (FEER) method to complex eigenvalue 
analysis in NASTRAN Is described. As In the case of real eigenvalue ana)’’sis, the elgen- 
solutlons closest to a selected point In the elgenspectrum are extracted from a r’educed, 
symmetric, tridiagonal elgenmatrlx whose order is much lower than that of the full-size 
problem. The reduction process is effected automatically, and thus avoids the arbitrary 
lumping of masses and other physical quantities at selected grid points. The statement 
of the algebraic eigenvalue problem admits mass, damping and stiffness matrices which 
are unrestricted in character, i. e. , they may be real, complex, symmetric or uns3nnmetric, 
singular or non-singular. 




rOREVVORD 


The Tridiagonal Reduction mothod for real eigenvalue analysis (Reference [l] ) 

is extended in tlie present report to accommodate complex eigenvalue problems; as 

t^Tpified by complex, unsymmetx'ic matrices. The implementation of this teclinique in 

NASTRAN pex’mits the extraction of specified numbers of eigenvalues lying closest to 

selected points in the complex plane, as well as the associated eigenvectors. The 

♦ 

eigenvalues may include multiplicities. 

Detailed descriptions of the theoretical development, computational procedures, 
and NASTRAN user instructions are provided. 
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INaU^ODUCTION 


Tlie complex Tridiagonal Reduction method Is im extension of the FEER nlgoritlim 
(Fast Eigenvalue Extraction Routine) for real eigenvalue analysis to complex, algebraic 
oigonproblem formulations. A specified number of eigenvalues lying closest to a selected 
point in the complex plane are sought, as well as the associated eigenvectors. As in the 
case of real eigenvalue analysis (Reference [1] ), these eigensolutlons are extracted from 
a symmeti’ic, tridiagonal eigonmatrix whoso order is much lower than that of the full-size 
problem. In fact, the size of this canonical, reduced mati'ix is of the same ox’der of mag- 
nitude as the number of desired I’oots, even if the discretized system model possesses 
tliousands of degrees of freedom. The reduction process is carried out via an automatic 
algorithm requiring a finite number of stops. Thus, a basic wealmess of methods requir- 
ing the lumping of masses and other physical quantities at arbitrarily selected degrees of 
freedom (Refei’cnces [2] - [4] ) is avoided in reducing the problem size. 

With regard to computational speed, the complex Tridiagonal Reduction method is 
somewhat slower tlian the Ilessenberg method (References [5] and [6] ) for small problems 
(on the order of one hundred or less degrees of freedom), if aR the existing eigensolutions 
are to be calculated. However, it becomes more efficient than the Hessenberg method 
when the number of requested eigensolutions is much less than the foil problem size. More- 
over, for much larger problems, the central memory requii’ement of the Hessenberg method 
exceeds the capabilities of most large computers, so that it becomes unavailable as a solu- 
tion option. This limitation does not exist in the case of the Tridiagonal Reduction Method. 

As shown in Section 2.1, tlie complex Tridiagonal Reduction method employs a 
’ single initial shift point, and hence only one matrix decomposition is required for each 
neigliborhood chosen in the complex plane. It therefore is more efficient than the complex 
Inverse Power method, which typically performs many shifts and decompositions for each 
region selected. 

The theory and computational procedux’es for complex analysis depart from those 
of real analysis in the following major respects: 

1. Both left and right bi-orthogonal vectors must be created in the process of 
constructing the reduced tridiagonal matrix. 
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2. The recliicod tridiagonal matrix, while symmetric in form, Is, in general, 
complex rather than real. 

3. The calculated theoretical errors in the computed eigenvalues are estimates 
rather than upper bounds. 

4. Eigensolutions closest to one or more specified points (shift points) in the 
complex plane are found. All eigensolutions obtained for previous shift 
points are swept out of the problem to prevent their regeneration when 
dealing with the current shift point. 

The development of the method is carried out in Section 2. A detailed summary 
of the computational procedures used in NASTRAN, flow diagrams and user Instructions 
are provided in Sections 3 and 4. NASTRAN user messages and optional diagnostic 
output are described in Section 5. 
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2 . 


THEORY FOR COMPLEX EIGENVALUE ANALYSIS 


2. 1 Problem Fonnulntion 

The gonernl complex eigenvalue pi’oblom is stated in the form 

[Mp^ +Bp •l■K]C^) = 0, (1) 

where [M] , [B], and Ck] may be X’eal, complex, symmetric or unsymmeti-ie, singular 
or non-singular. A specified number of eigenvalues, p, lying closest to a specified point, 
(called a sliift point) in the complex plane are to be found, as well as tlie associated 
eigenvectors fiij . The eigenvalues may include multiplicities. 

A velocity vector, fv} , Is defined as 

{v}=pCu3, (2) 


and a shifted eigenvalue, X, Is defined as 


x=p-x„. 


( 3 ) 


Substituting Equations (2) and (3) into Equation (1), the result may be written in partitioned 
form as 



* 
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or, using the inverse form, results in 


where 


[A][x} = ACx} , (5) 
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( 7 ) 


nnd 



( 8 ) 


Equation (5) is an eigenvalue problem in standard form, showing that the order of 
the eigenvalue problem is doubled due to the presence of the [B] matrix. In the special 
case where [b] is null (e. g. , no damping) , the problem formulation becomes 

[Mp^+K] Cu} =0, (9) 

and the double-size eigenvalue problem can be avoided by considering the mathematical 
2 

eigenvalue to bo p . Thus, let 


and 



( 10 ) 


( 11 ) 


Substituting Equations (10) and (11) into Equation (9) and using the inverse form results in 


[K +x\l] (u] = A {u} . (12) 

Comparing the above with Equation (5) shows that the standard form with a null [b] matrix 
is acliieved by setting 

[A]=Ck+X^M]“^[-M], (13) 

and 
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Cx} « Cu} . 


(14) 


> 


Stnco tho elgenmatrlx, [A}, Is, In goncvnl, imsymmotHc, the olgenvcctors, (xj, are 
{jrthogonal to tho eigenvectors, (x}, of tlio transpose eigonproblem 

[A]'^Cx] » A{x3 , (15) 


so that for A. A, , 

1 j 

Cx.}^tXj3 = 0; (1C) 

The above relationship is a biorthogonality condition and the associated eigenvectors, (x.) 

and {x, 3 , are called right and left eigenvectox’S, respectively. 

J 

2. 2 The Reduction Algorithm 

A reduction of the order of the eigenvalue problem, Equation (5), is effected through 
the transformation 


{x3 = [V] {y3 , (17a) 

nxl nxm mxl 

and 

Cx3 = [V] {y} , (17b) 

nxl nxm mxl 


where {x} and {x3 are approximations of (x) and Cx}, respectively, n is the order of 
die imreduced problem, and m s n. The above transformation matrices are chosen to be bi- 
orthonormal, so that 


[v]'^[v] = m. 


( 18 ) 


From JHquations (5), (17), and (18), it is seen that 
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whore 


ClI] Cy3 = Afy] , 


(19) 


Di3'=I:v]'’^[a] [V], (20) 

mxm 

and A- is an approximation of the elgenvaluo, A • 

Thus, Equation (19) Is an m^*^ order eigenvalue problem, where ms:n. The value 
of m is established according to the criteria given in Section 2.3. 

As in the case of real eigenvalue analysis (Reference [1] ), the Lancv.os algoi'ithm 
is used to construct the transformation matrices vector by vector, 1. e. , 

Cv] « [{v^} , Cvg) 

nxm 

[V] = C{v^3 , Cvg} ,•••♦. C v^3 ] , (21b) 

nxm 

such that the reduced mxm matiix, [ll] , is tridiagonal and its eigenvalues accurately 
appx'oximato the roots of Eqviation (5) having the largest magnitude (or, equivalently, the 
physical roots, p, closest to the specified point of interest, X^, in the complex plane). 

The form of the algoritlmi for generating successive vector paii-s is, according to 
the Lanezos technique, 


,T,= , = 



(22a) 

■"i, 2 ' V- • 

(22b) 


where and (v^3 are starting vectors (see Section 2.4), l:sism-l, and all m 

vectors pairs are biorUionormal. TMs implies that 

a ={v3^[A]fvJ, 

If J J ^ 

-G- 


(23a) 



”l,j ° » 


while the normuli'Mng faelors are given by 




Using symmetry arguments similar to those employed for real eigenvalue analysis 
(Reforenee [1] ), it ean be shown that 


“i»i ’ 






The tr^iVifr »med, reduced eigemnatrix in Equation (19) is consequently tridiagonal 
and symmetric, having the form, 


■W = [V]'-^CA] [V] = 


“22 dg 

■■s “>33 ‘’4 

N \ 

S 'v 

V N M 


S V 

''d ^ a 

m mm 


The matrix coefficients are theoretically given by the simplified recurrence formulas 
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^[aKVj) 


'Vi> ° '''i’ ■ ‘'i.i'''i^ ■ 


L 1^1, m 




(27a) 




y i = 1, m-1 (27b) 


where the sequence is initialized by choosing random, biorthonormal starting vectors for 
(Vi), [v^J and by setting d^ = 0; {v^j = {v^ } = Co) . 

The final off-diagonal term, d _ , given by Equation (27a) is used in establishing 

Itl+X 

error estimates for tlie computed eigenvalues (see Section 2. 6), In addition, the above 
algorithm is modified in the computational scheme as follows: 

1. Each pail* of vectors [v,j^_3 > calculated in Equations (27b) is reor- 

i*rl 1 t"X 

tliogonalized to all previously computed pairs, as described in Section 2. 5, 
before re-entering Equations (27a). 

2. The size, m, of the reduced problem is a function of the number of accurate 
eigenvalues requested by the user and is limited to the number of ftr Ue 
physical eigenvalues available (see Section 2.3). 

The eigenvalues. A, and eigenvectors, {y} , of Equation (19) are extracted using the 
Q-R iteration algorithm and eigenvector computational scheme described in connection with 
the Upper Ilesscnberg method in NASTRAN (Reference [5] ). They are then converted to 
physical form as follows: 
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[V](y^ 


y [B] ■/ Co] , (28a) 


tuj) = Cv] Cy,) 


^ [B] = Co] . 


(28b) 


The velocity vector Cv^,} in Equation (28a) is discarded prior to further processing of the 
set of eigensolutions by NASTJRAN. Also, all solutions which fail the FEER error test are 
rejected. However, the number of acceptable solutions will, in all probability, equal or 
exceed the number requested by the user when the reduced problem siae is chosen accord- 
ing to the criteria described in the following section. 

2, 3 Criteria for the Size of the Reduced Eigenvalue Problem 

The maximum number of finite eigensolutions, including any existing rigid body 
modes. Is equal to th . rark, r, of the elgennmtrlx, [A], In Equation (5). Thus, for example, 
massless degrees of freedom, appearing as zero diagonal terms in the [M] matrix, will 
. result in singularities (ranlc reduction), which imply infinite physical eigenvalues. These 
spurious roots are swept out of the problem in the complex FEER process (see Section 
2.4) with a consequent reduction in the available eigensolutions. 

A further consideration in limiting the maximum problem size is that the user has 

the option of requesting eigensolutions in the neighborhood of several shift points 

(X X , . . . ) in the complex plane. In the Tridiagonal Reduction method, all eigen- 
01 0 ^ 

solutions, f, obtained for previous shift points are swept out of the problem to prevent 
their re-generation when dealing with the current shift point (see Section 2. 5). This implies 
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that the maximum possible sIko, m, of the reduced problem is further limited to 


m 

max 


X’ - f . 


( 29 ) 


On the basis of numerical experiments, similar to those cited in Reference [1] for 
real eigenvalue analysis, it has been found that when m«m a first grouping of more 
than m/2 computed eigenvalues closest to the shift point are in accurate agreement with 
the corresponding number of exact eigenvalues, provided that 7 Sm ^ ^^max’ I’emain- 
ing reduced-system roots ax'o spread across the remaining exact eigenspectrum. To en- 
hance the accuracy of the associated eigenvectox’S, the minimum problem size is further 
increased to twelve, again assuming that m <s: m 

max 

Thus, if the user x’equests a total of q eigenvalues closest to a specified point in 
the complex plane, the order of the reduced problem is initially set to 


= min [(2q+10), (2n-f)] ; 

Cb] * [0] . 

(30a) 

= min [ (2q+10) , (n-f)] ; 


(30b) 


Although the total number of eigensolutions requested should not exceed 

tliero is usually no simple way to discern this upper limit in complex eigenvalue problems. 

However, the reorthogonalization tests of Section 2. 5 are designed to automatically 

establish this upper limit. If the latter tests fail for some vector pair {v, [v. _ ) i 

1+1 1+1 

this is an indication that a null vector has been generated because m linearly inde- 
max 

pendent vectors have already been obtained. The recui’rence algoiatlim, Equations (27), 
is then temiinated and the order of the eigenproblem is further reduced to m = i. 

2.4 Choice of the Initial Ti’ial Vector and Restart Vectors 

Because of the inverse relationship between the computed eigenvalues, A, and the 
physical eigenvalues, p, (see Equations (28)), spurious eigenvalues cori'esponding to 
A = 0 are equivalent to p . Since these eigenvalues anc. their corresponding eigen- 
vectors are of no intei’est and may cause nmnexlcal instabilities, they are eliminated from 
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the voducecl trlcUngonal problem by employing a constrained set of (v], Cv) vector pairs 
having zero projection on the set of eigenvectors associated with A~ 0. However, any 
non-null right and left vectors {w^^l , [Wj^j , genoratod from any other non-null vectors 
{Wj.} » Cw^,3 thx'ough the transformations 

{wJ=[A]{w^3. (31a) 

(w^}« [A]*^Cw^} , (31b) 

will contain no components of the eigenvectors corresponding to A = 0. Further, after the 
above vectors are orthonormalized 




(32a) 


{v^} ^ rp *1 /o 3 } 


(32b) 


the next vectors generated, 

(Wg) = Ca] Cvj^3 - 

(Wj) = CA]'^t?j_) - . (33b) 

as well as all subsequent trial vectors will be free of A = 0 eigenvectors. 

Thus, the starting vectors, {v^3 » Cv^3 , in the reduction algorithm are generated 

from Equations (31) and (32), starting with pseudo-random vectors {w^3 and (w^,3 . If 

these pseudo-random vectors should, by chance, be deficient in some true eigenvector 

components, then a null vector W. ^'^i+l^ may be generated at some point in 

Equations (27). Tliis is indicated by the appearance of an off-diagonal term, d.^^, wMch 

is exceedingly small compared to the corresponding diagonal term, a. in the reduced 

1 # 1 
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tvicllagonal matrix. Tho tost used is that a null vector has been generated when 





( 34 ) 


where t is the number of decimal digits carried in tho computations. In this event, new 

restart vectors are employed for . These vectors are generated exactly 

as in the case of tho initial tx'ial vectors, but using difforeait pseudo-random number seeds. 

The recurrence algorithm for generation of subsequent vector pairs (v. ^3 » [v ,3i etc., 

iHT iHT 

is then continued in the usual manner tmtil the required number of vector pairs has been 
generated. 


2. 5 Sweeping-out of Previously Obtained Eigenvectors and Reortliogonalization of the 
Ti'ial Vectors 


As in the case of real eigenvalue analysis, successive trial vectors generated by 
the theoretical algoritlim. Equations (27), tend to degrade rapidly as the computations pro- 
ceed in the finite digit calculations, such that later right vectors are far removed froxn Ortho- 
gonality to earlier left vectors, and previously computed eigenvectors, {x ] , (x 3 ? • <- r ^ , 

— _ X u 

(xj; {x 3» {x„3 »••••> Cx 3 * Therefore, each new vector pair fv. ,3 > [v. ,3 obtained 

1 X ^ I X+1 iH'l 

from either the algorithm or the pseudo-random number generator (see previous Section) 
is denoted as {v.f^|3 > and reox’thogonalized with respect to all the previously ob- 

tained vectors. Tliis is accomplished by iterating, 


cv - A -£ ctijf txj] . 


i+1 


3=1 


3=1 


3 3 1 


(35a) 


(35b) 


where k = 1, 2, . . . , s mitil an acceptable vector pair [v|®^3» found which satis- 

fies the orthogonality criteria 


max 
1^3^ i 




(36a) 
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max 




max 




max 



(5GjjJ 

(3Gc) 

(3Cd) 


where s s lo, and t is the number of digits cai’ricd in the computations. If, for some 
value of i, the above criteria are not satisfied, then it is assumed that a new pair of 
trial vectors cannot be generated and a reduced eigenvalue problem of order m = i is 
solved as discussed in Section 2. 3. 

If the above criteria are met, then the resulting vectors are biorthonormali'/ied and 
set equal to the new trial vector pair, l.e. , 


fv 1= 


(37a) 


Cv,^J = 


1 . 


tv '!>) . 




(37b) 


These new vectors are used to compute the next off-diagonal tenn in the reduced tridiagonal 
matrix from the formula 


I^Vi) 


However, if 


(38) 




(39) 


it is probable that either or is a null vector. Di this event, a reduced 

eigenvalue problem of order m = i is solved, as above. 
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If the criteria given by Equations (30) and (39) are both 3uet, then the new, norma- 
lized vectox’s are used to continue the reduction algorithm. 

2.0 Error Estimates for the Computed Eigenvalues 

Followhig a developmant similar to that of Reference [l] for real eigenvalue 
analysis, It can be shown that 

l|Ai|-|Aj|s|Viy.J- (AO) 

The above show« that the absolute value of the difference betiveen the computed and tnie 
eigenvalue magnitudes is propox’tional to the magnitude of (wliich is the next off- 

diagonal term that would be generated had the reduced tridiagonal mati’ix, [H], been in- 
cx’eased from order m to order m+1) and (which is the last term in the reduced- 
system eigenvector associated with A.)« 

Converting Equation (40) to physical eigenvalue form, using Equations (8), (10), 
and (11), yields, 



(41a) 


(41b) 


The use of the above erx'or bomid estimates as critex’ia for acceptable eigensolutlons 
is desexibed in Section 3, 
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3. 


SUMMARY OF COMPUTATIONAL PRQCBDUIIES AND FIX)W CHARTS 


Flow diagram? illustrating the computational procedures are shown in Figures T 
and 2. The details of each block are summarised below. 

(1) Establish Tentative Reduced Problem Slso (see Section 2. 3) 

Calculate a tentative size, m, of the reduced eigenproblem from 


. =min[(2q+10),(2n-f)]; [B] 7 ^ [O] 

^ =min [(2q+10), (n-f)]; [b] = [o] , 


(la) 

(lb) 


where n is the ox’der of the [K] > [B], and [m] matrices, q is the number of user re- 
quested eigensolutions , and f is the total number of previously extracted eigensolutions 

2 

calculated for earlier shift points (X or X ). 

0 o' 

(2) Construct Factors of Dynamic Matrix 
(a) Set 


p = [K +X +X ^ M] ; [B] ^ [O] 

Cd] 

^ ; [B] = C0]. 


(2a) 


(2b) 


(b) Deco3npose the above dynamic matrix: 


Cd]-[L][U], (3) 

where Cl] Is a unit lower triangular factor and CuJ is an upper triangular far'or, 
using complex arithmetic (subroutine CDC0MP in NASTRAN). 

Save the triangular factoi’s for later use in premultiplication opei’ations involving the 
eigemnatrix, [A] , and its tx’anspose. If the decomposition fails, the sliift value is 
increased by A = .02(l+i) and a new decomposition is attempted. If tliis second 



clecoinposUton is unsuccessful, then It is assumed that the problem contains unremov- 
able slngularttlcs In the nejgliborhoocl of the current slilft point and the next I’cglon, if any, 
Is examined. 

(3) Execute Complex Tridiagonal Reduction Algorithm (see flow diagram for this block, 
Figure 2) 

3.1 Initialise the Ilocurronce Algorithm (see Section 2.2) 

Initialize the vector index to i=0 and set 


(v^) = (? ) = to) . 

o o 


(4) 


3. 2 Generate a Pair of Starting or Restart Vectors and set d. = 0. 0 (see Section 2.4) 

l*rX 

a. Construct right and left vectors fw^] and {w^3 using a pseudo-random niunber 
generator. The vectors contain 2n elements for [B]?^[0] and n elements for 

Cb]=[o]. 

b. Sweep-out components corresponding to A = Oj 


f'Vi3 = . 


t'Vl) = 


(5n) 

(5b) 


c. Normalize the above vectors; 


(v<")) = 

’• i+l' 


T. ,4/2 tVl’- 


Ct'VHx] Cw^,,)J 


(6a) 


7- fw 1 

f'Vi” 


[?W) = 

i+1-' 


(Gb) 


d. Set d^^^ = 0.0 and proceed to block 3. 5. 
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3, 3 Crcnio Ono Pnir of Approximate Trial Vectors nncl One Dingonal Coefftclont 
(see Section 2.2) 

The rcciiri’oncc algorithm is: 


ni^i=Cvjf[A] Cvp, 

CWj^l3 “ UKv^)- a.^^(v^3-cljCv._^3 , 
CWi.^l3 = [A] c Vj3 - jC V.3 -• d^C Yi3 > 




“l+l 






d-xi 
1+1 

where the above vectors are approximations to the new right and left trial vectors. 

3. 4 First Normalization Test (see Section 2. 4) 

The test is: 


( 7 ) 




( 8 ) 


Pass: Proceed directly to block 3. 5. 


Fail: Return to block 3. 2, generate a new pair of restart vectoi’s for Cv/^^3 and 

= tOl 

^ i+P’ and proceed to block 5.5. 

3. 5 I terate to Obtain OrthORonalized Trial Vectors (see Section 2. 5) 

Designate {x,3, [x.3 (j = l» as previously calculated and stored eigenvector pairs, 
1. e. , calculated eai’lier by complex FEER for previous shift points in the complex plane. 
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Perform the Ucrntlons, 




j=l ^ ^ 


(fla) 




J=l 


y ^ m 


k - 1,2. , . , until 


- s fi,! i 


j=l 


y "1+1 


l[vfCv^''hUlO^"‘^ 

i^j^i '■Vi^' » 

max I c ,T(.= (k), ) 2-t 

uj^iUV ' 

max I f= ,T. {k)-, I . . ,2-t 


max tr (^)i i 




m 


(10) 


or k - 10. 

If the orthogonality criteria , Equations (10), are satisfied, proceed to block 3.6. Other- 
wise, set the problem size, m, equal to i and proceed to Exit, 

3. G Normalize the Orthogonalized Trial Vectors (see Section 2. 5) 

Compute 




tv'+J,}. 


(lla) 
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These are the new orthogonali/^cd and norinulisicd trial vectors. 


3 . 7 Second Norinalization Tost and Ci’oation of OCf-DIngonal Coefficient (sop Section 2,5) 

(a) Compute the next off-diagonal term of the reduced tridiagonal matrix from 

(b) Verify whether the following test is mot: 

(I®’ 

If it is, set i = i+1 and return to block 3. 3 for continuation of the recurrence 
algoi-ithm. If the test falls, set m=i to reduo© the problem si»e and proceed 
to Exit. 

Solve Reduced-System Eigennroblem (sec Section 2. 2) 

(a) The coefficients a, , , a , . , . , a and d , d d computed in block 

XX 22 Itllll 2 o 111 

3 form the diagonal and off-diagonal terms of an m ox'der symmetric, tri- 
diagonal matrix, Ch], (the matrix is, however, usually complex rather than real), 
til 

(b) The m order eigenvalue problem 

Ci-I] Cy} = A ( y) (14) 

is solved for the eigenvalues. A., and the eigenvectors, CV-l » using the Q-R 
iteration algoritlun and eigenvector computational scheme described in 
connection with the Upper Ilessenberg method in NASTRAN (Reference [5]). 

(c) The reduced system eigenvectors are normalized so that 
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( 15 ) 




(5) Computo EsUmalo of ElKonv»tue Errors (see Section 2, 0) 


(a) Estimates of absolute relative errors In the computed eigenvalues are obtained 
from 




t3 


IfilM 


1 


« 


lA.I 

& 


;[b]5^Co3, 


(IGa) 
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[B] = [0] , 


(1Gb) 


where d^^^^ Ib the last off-diagonal terni computed in block 3. 3 and Is 
the last clement in the vector tyj^} . If the physical eigenvalue, p^, corresponds 
to a zero root (e.g. , a rigid body mode), the above computational scheme is 
invalid and therefore bypassed. A zero root is assumed to occur whenever 


h‘il 

Kills 



I 


where, 

KMs = i + + 


(17) 


(18) 


mid is denoted by setting the error to zero, 

(li) The eigenvalues are listed in order of increasing distance from the shift point, 

X , to determine whether their associated estimated errors, meet an 
o 1 

acceptable relative error tolerance set by the user on the EIGC bulk data card 
(the default value Is O.lO/n, where n Is the order of the stiffness matrix). 

The first eigenvalue not meeting the tolerance test, as well as all subsequent 
eigenvalues furtlier removed from the center of interest, are considered to lack 
sufficient accuracy and are therefore discarded. 
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(c) Acceplnblo oigonvuluos obtainca In the above manner are reordered accordijig 
to the magnitude of the Imaginary part, with positive values considered as a 
group ahead of all negative values. 

(Q) Compute Physical Eigenvalues and ElKonvectors and Store (sec Section 2. 2) 

. Tho mathematical eigenvalues, A^, and eigenvectors, are converted to 

physical form as follows; 




V Co] 


"'1 =Cv]Cy^3, 

i 


(19a) 


where 


p. lm(Pj^)>0 

' ^i 
Cu.3 = CV] {y^l 


r [B] = Co] 


CvD^CCvj^). Cvg), ,CvJ]. 


(19b) 


(20) 


The velocity'vector [v^^ } in Equation (19a) is discarded. 
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Figure Overall Flow Diagram - Complex Tridiagonal Reduction Method 
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Figure 2, Flow Diagram for Block 3, Execute Complex Tridiagonal Reduction Algoritlun. 
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NASTRAN USER'S INSTRUCTIONS 


Tho following pages show modifications of the EIGC card in the NASTRAN bulk 
data deck which accommodate user implementation of the Tridiagonal Reduction method 
for complex eigenvalue analysis. The modifications consist of additions to the standard 
user instructions and are underscored for ease In Identification. 

When the complex Tridiagonal Reduction method is invoked, the E parameter on 
this card represents the maximum allowable value of the computed absolute relative 
error in a physical eigenvalue (see Section 2.G). If this value is exceeded, the associated 
eigensolution is not accepted for further processing by NASTRAN. A detailed list of the 
maximum relative errors computed by complex FEER can be obtained by requesting 
DIAG 12 In the NASTRAN Executive Control Deck (see Section 5). 
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BULK DATA DECK 

Input Data Card EIGC Complex: Eigenvalue Extraction Data 
Description: Defines data needed to perform complex eigenvalue analysis 


w 



Format and Example: 


, 1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

EIGC 

SID 

METI-l(2fD 

NC?JriM 

G 

C 

E 



+abc 1 

EIGC 

. 14 

DET 

PdlNT 

27 


1.-8 


* 

ABC 


+abc 

ttal 



ll 

•*^el 

msEsm 


-t-def 

+BC 

HOQIIII 

5.6 2.0 


2.0 

4 

4 

1 



+def 


mssBm 

“b2 

^^b2 

^2 

— 

HEM 


+EF 

Eu 

-5.5 

5.G 

5.6 

1.5 

6 

3 

^ ^ ! 


(etc. ) 
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Field Contents 

Sot identification number (unique integer > 0) 

Method of complex eigenvalue extraction, one of the BCD values, "INV”, 
"DET", "HESS" or " FEER " 

INV “ Inverse power method 

DET - Determinant method 

HESS “ Upper Hessenberg method 

FEER - Tridiagonal Reduction Method 

Nt^RM Method for nonnalizing eigenvectors, one of the BCD values "MAX" or 

"P(2JlNT" 

MAX - Normalize to a unit value for the real part and a zero value 
for the imaginary part, the component having the largest 
magnitude 

PC^INT - Normalize to a unit value for the real part and a zero value 
for the imaginarj' part the component defined in fields 5 and 
6 “ defaults to "MAX" if the magnitude of the defined com- 
ponent is zero. 

G Grid or scalar point identification number (Required if and only if 

N0RM=P0INT) (Integer > 0) 

C Component number (Required if and only if N0RM="PC!lNT" and G is a 

geometric grid point) (0 ^integer s 6) 


SID 

METIhJ^D 


E Convergence criterion (optional) (Real ^0. 0) 

For method = "FEER", error-tolerance on acceptable eigenvalues in per- 
cent (default value is .lO/n, where n is the order of the stiffness matrix) 


<%■ ] 

<“bj’ “bj> ] 


) 


Two complex points defining a line in the complex plane (Real) 

For method = "FEER", is a point of intei-est in the complex 

J .1 - ‘ i • t .*1 - ] 1 1 ..!. 111 .. 1 / 


plane, closest to which the eigenvalues are computed; 
The point w,^j) is ignored. 


a 


2 ] 


+ to 


aj 


0 . 


Widtli of region in complex plane (Real >0.0) 
Blank for method = "FEER". 
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Estimated numbex’ of roots in each region (Integer > 0) 
Ignored for mctliod ~ ’*FEER" . 




Desired number of roots in each region (Default is 3Ngj) (Integer >0) 
Desired number of accurate roots for method ~ "PEER'' (Default is 1) . 


Remarks ; 


1. Each continuation cax’d defines a rectangular search x’egion. For method = 
"PEER", tlie card defines a cii'cular search region, centered at (o;nj» WajL 
and of sufficient radius to encompass N^j j roots. Any number of regions 
may be used and they may overlap. Roots in overlapping regions will not be 
extracted more than once. 

2. Complex eigenvalue extx’action data sets must be selected in the Case Control 
Deck (CMETH0D=SID) to be used by NASTRAN. 

3. The units of O', to and -6 are radians per unit time. 

4. At least one corfinuation card is required. 

5. For the dp' -rminant method with no damping matrix, complex conjugates of 
the roots found are not printed. 

6. See Section 10,4.4. 5 of the Theoretical Manual for a discussion of convergence 
cx'iteria. 

7. For the Upper Hessenberg metliod, controls the number of vectors com- 

puted. Only one continuation card is considered and the (a, to) pairs, along 
with the parameters and Nqi, are ignored. Insufficient stox*age for HESS 
will cause the program to switch to INV. 

8. The eri'or tolerance, E, for the "FEER" method is with regard to 



for Cb] [0] and 


for [B] = [0] , 


where p^^ is a computed eigenvalue and p. an exact eigenvalue. 


- 27 > 





V 


5. USER MESSAGES AND OPTIONAL DIAGNOSTICS 

5. 1 NASTRAN Function Module User Mossages for the Complex Tricliagonal Reduction 
Method 

5.1.1 General 

The following is a description of the NASTRAN user messages wliich may be 
generated bj' NASTRAN during the execution of the Complex Tridiagonal Reduction method 
and which are unique to this method. Explanatory information is pi’ovided following the 
text of each message and, in the case of a fatal message, corroctiv6 action is indicated. 
Refer to the NASTRAN Users' Manual, Section G for a complete listing of other system 
and user messages. 

Fatal messages cause the termination of the execution following the printing of the 
message, text. These messages will always appear at the end of the NASTRAN output. 
Wai’ning and information messages will appear at various places in the output stream. 

Such messages convey only warnings or infoi’mation to the user. Consequently, the exe- 
cution continues in a normal manner following the printing of the message text. 

5. 1. 2 List of User Messages 

3149 *■** USER WARNING MESSAGE 3149, USER SPECIFIED NEIGHB0RH(2i0D 

CENTERED AT (2^RIGIN N0T ALL(2JWED, CENTER SHIFTED T(^ 

THE RIGHT . 001. 

Point of interest in the complex plane closest to wMch the 

eigenvalues will be computed, was input as (O.'o, 0.0) on an EIGC bulk 
data continuation card. Since tins is an inadmissible choice, the point 
automatically used was (.001, 0.0). 

3150 ■*** USER WARNING MESSAGE 3150, DESIRED NUMBER 0F EIGENVALUES 

**** INVALID. SET = 1. 

Number of accuj*ate roots desired was omitted, input as zero or 

negative on an EIGC bulk data continuation card. The number auto- 
matically used was 1. 

3151 USER WARNING MESSAGE 3151, DYNAMIC MATRIX IS SINGULAR 

(Occurrence **=•'*) in neighborhood centered at 

Point of interest in the complex plane (oijjj^, closest to wliich the 
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3152 *** 


3153 


3i54 


3155 *** 


3156 *** 


eigenvalues will be computed, was input (oo close to an eigenvalue 
on an EIGC bulk data continuation eax'd, Tlie point is automatically 
shifted by adding , 02 to both the real and imaginary parts. If llie 
dynamic matrix is still singular, the next neighborhood, if any, is 
searched. 

USER INFC^RMATIC^N MESSAGE 3152, SUBRC^UTINE ALLMAT (J^UTPUT 
EIGENVALUE **** IS NULL. 

When an eigenvalue output from subroutine ALLMAT is exactly Sicro, 
the formula for computing the associated theoretical error test fails. 

The magnitude of the eigenvalue is considered to be 10”^® for use in 
that formula. 

USER WARNING MESSAGE 3153, ATTEMPT T0 NORMALIZE NULL 
VECT(2iR IN SUBRC^UTINE CFEER4. N0 ACTI0N TAKEN. 

An eigenve« ^r output from subroutine ALLMAT is a zei^o-vector. 

USER WARNING MESSAGE 3154, SIZE 0F REDUCED PRQiBLEM 
DECREMENTED 0NCE (NQ^W =«=*+*) DUE T0 NULL ERR0R ELEMENT. 

If subroutine CFEER4 receives a reduced tridiagonal matrix having 
error element dj^^.j^ exactly (0,0), it is impossible to compute mean- 
ingful theoretical error estimates for any of the eigenvalues. The size 
of the reduced problem is reduced by one, so that d^^ becomes the new 
error element. 

USER WARNING MESSAGE 3155, REDUCED PROBLEM HAS VANISHED. 
N0 R00TS F0UND. 

If decrementing the size of the reduced problem (see message 3154) 
causes the size to become zero, the program continues to the next 
neighborhood, if any. 

USER WARNING MESSAGE 3156, SIZE 0F REDUCED PROBLEM 
RESTORED T0 ***+ because NEXT EER0R ELEMENT WAS ALS0 
NULL. ERR0R ELEMENT SET = **** 

This message follows message 3154. If dj^ is also exactly zex’O (in 
addition to being exactly zero), then the ox’iginal reduced px’oblem 

size is restored and dj^.j .2 is set to (£, 0) whex'e & == E/lOO and E is 
the error tolerance on acceptable eigenvalues input on the EIGC bulk 
data card. 


ORIGINAL PAGE iS 
-29- OF POOR QUALITY 



4 


3157 *** 


3158 *** 


3159 *** 


3160 *** 


3161 


3162 **>(! 


3163 **♦ 


USER WARNING MESSAGE 3157, PEER PR^^JCESS MAY HAVE 
CALCULATED FEWER ACCURATE M(^DEfI **** THAN REQUESTED 
IN THE NEIGHB^RIK^CJD C^F **** **** 

Tlio desired number of eigenvalues specified on the EIGC bulk data 
continuation card exceeds the additional number that can be calculated 
by the Complex Tridiagonal Reduction (Complex FEER) method in the 
current noighborliood, 

USER WARNING MESSAGE 3158, NCiJ ADDITIO^NAL M0DES CAN BE 
F0UND BY FEER IN THE NEIGHBCiRII00D QF 

An initial pseudo-random vector cannot bo made oi’thogonal to the exist- 
ing set of orthogonal vectors (which come from Restart and from all 
prior-neighborhood sets of eigensolutions). 

USERINFO^RMATIQN MESSAGE 3159, ALL S6<LUTI(2iNS HAVE BEEN 

fQund. 

The FEER method has solved the entire problem. Any additional neigh- 
borhoods (as specified by the presence of EIGC bulls data continuation 
cards) are ignored. 

USER INFOJRMATI0N MESSAGE 3160, MINIMUM QPEN CC^RE N(?5t 
USED BY FEER **** W0RDS {**** K BYTES). 

This message indicates tlie amoimt of open core, in both b 3 rtes and 
words, not used by FEER. 

USER WARNING MESSAGE 3161, DESIRED NUMBER 0F EIGENS0LU- 
TldNS FCIR NEIGHB0RHCj(2^D O^F ***'* CENTERED AT 
**** **** EXCEEDS THE EXISTING NUMBER ****, ALL EIGENSC2JLU- 
TK2INS WILL BE SOUGHT. 

The desired number of eigenvalues specified on the EIGC bulk data 
continuation card exceeds tlie size of the eigenmatrix, which is the 
maximum possible number of existing eigenvalues. 

USER WARNING MESSAGE 3162, ATTEMPT TC^ N(2lRMALIZE NULL 
VECTOR. NQi ACTICiN TAKEN. 

The general vector normalization routine (CFN0R1 or CFNQR2) has 
a zero-vector input to it. 

USER WARNING MESSAGE 3163, ALL **** SC)LUTI0NS HAVE FAILED 
ACCURACY TEST. NQ5 RC^^TS FQUND. • 
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Tlie number of oigcnsolutions passing the I’elativo error tost is sioro, 
The maximum allowable error Cor the relative error lost is specified 
in field 7 of the EIGC bulk data card. A detailed list of t)ie computed 
' error bounds eould have been obtained by requesting DIAG 12 in the 

Executive Control Deck. 

31C4 USER INF(2^RMATI0^N MESSAGE 31G4, ALL >»=*** Sf^LUTIONS ARE 

ACCEPTABLE. 

Ail the eigensolutions obtained in the reduced problem corresponding to 
the point of interest pass the relative ox’ror test. The maximum ailow- 
able error for the relative error test is specified in field 7 of the 
EIGC bulk data card. A detailed list of the computed error estimates 
could have been obtained by requesting DIAG 12 in the Executive 
Control Deck, 

31GG *** USER INFC2lRMATIC2lN MESSAGE 3165, **** SCSLUTICJNS HAVE BEEN 
ACCEPTED AND **** S0LUTI0NS HAVE BEEN REJECTED. 

Some eigensolutions passed the relative error test and some did not. 

31G6 *** USER INF(^RMATI0N MESSAGE 31G6, ****■ MC^RE ACCURATE EIGEN- 
SCJLUTI0NS THAN THE **** REQUESTED HAVE BEEN FQUND FQR 
NEIGHB(2iRII(2f<2iD **** 0F **** CENTERED AT **** ***=1^ USE DIAG 
12 TOJ DETERMINE ERR0R ESTIMATES. 

The number of eigensolutions passing the relative error test is greater 
than the number requested on the corresponding EIGC buUc data con- 
tinuation card. The maximum allowable error for the relative error 
test is specified in field 7 of the EIGC bulk data card. A detailed 
list of the computed error estimates could have been obtained by re- 
questing DIAG 12 in the Executive Conti’ol Deck. 
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5,2 Tho Eigenvalue Summary Table 

The following summary of tho eigenvalue amilysls performed, using the complex 
Tridiagonal Reduction (FEER) method, is automatically printed; 

1. Number of eigenvalues extracted. 

2. Number of starting points used. 

This corresponds to the total number of random starting and restart vectors used 
by the complex PEER process for all neighborhoods. 

3. Number of starting point moves. 

Not used in PEER (set equal to zero). 

4. Number of tr-iangular decompositions. 

Always equal to tiie number of points of interest (neighborhoods) in the complex 
plane processed by PEER, since ordinarily only one triangular decomposition is 
required by PEER for each point of interest, imless the dynamic matrLx is singular 
at a given point of interest, in which case an additional decomposition is required 
(obtained by moving the point of interest slightly). 

5. Total number of vector iterations. 

The total number of reorthogonalizations of all the trial vectors employed. 

G, Reason for termination. 

(0) All, or more solutions than the number I’equested by the user, have been 
determined (normal termination). 

(1) All neighborhoods have been processed, but PEER has not obtained the de- 
sired munber of roots in each neighborhood, possibly because tliey have 
already been found in other neighborhoods. 

(2) Abnormal termination - either no roots found or none pass the PEER error 
test. 
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5. 3 Optional Diagnostic Output 


Tho user can obtain special detailed Iriformation relating to the generation of the 
reduced problem size, tho elements of tho reduced trldlagonnl matrix, vector rcortho- 
gonallzatlon itei’atlons, computed error estimates, order of eigenvalue extraction, and 
distance of oxtraolcd eigenvalue from tlie center of interest by requesting DIAG 12 In the 
NASTRAN executive control deck. 

The meaning of this iiiformatlon is explained below in the ox’dor in which it appears 
in tho DIAG 12 output. 

(FAST EIGENVALUE EXTRACTIQ^N RCJUTINE)**** 

Tliis header is always printed first, 

****SINGLE PRECISICiN W0RDS C^F 0PEN C(2 JrE NC5t USED (SUBRiiiUTINE XXXX) 

- Open core not used by subroutine XXXX, in eiugle-preclsion words. 

XXX;'! ~ Either CFCNTL, CFEER3, or CFEER4. This message appears three times. 

Cl O' .CL ACCURACY CRITERI0N-* (INPUT VALUE **) 

* - Accuracy criterion, used for rejecting eigensolutions (expressed as a percentage). 

4 k - Valuo of nccux’acy criterion input by the user on tlie EIGC bulk data card. 

CFCNTL NEIGIIB0RII0(2(D * CENTER = ** ** N0.DES.RTS.= *** N0NSYM = 

+ - Positive integer indicating which neighborhood, or center of interest, is currently 

being processed. 

** - Center of intex'est in the complex plane. 

- Number of desired roots for the current neighborhood, input by the user on the 
corresponding EIGC bulk data continuation card. 

- Indicator which, when nonzex*o, forces the program to consider the matrices as 
non-symmetric, even though they may actually be symmetric. Tliis is input by 
the user in field 7 of each EIGC bulk data continuation caxxl. Tliis input was 
used during program checkout of the complex FEER process, and it should have 
no affect on the solution. However, the user should leave field 7 blanlc on each 
EIGC continuation card. 
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RECillTHOiGO^NALT/iATIClN ITERATiCiN ♦ TARGET VALUE *= 

ERRC^RS *= *** *** *** 

+ " The reorlhogonalizatlon lleration-nmTiboi’. This nrossago will appejir many 

times, as the FEER proeoss "cleans up" each trial vector by forcing It to be 
as orthogonal ns possible to the set of vectors already computed, 

** - ‘ Convergence toleranee, such tluit the errors must bo smaller than this value. 

In oi’der to avoid taking square roots, the tolerance and errors are all squared. 

*♦* .. Four reoi'thogonalisatlon errors, the first two of which correspond to 
the orthogonality of the current right and left handed trial vectors, respectively, 
with respect to all previously computed vectors in the current neighborhood, 
and the latter two of which correspond to orthogonality of the same vectors 
with respect to all eigenvectoi’S previously computed (restart and prior neigh- 
borhoods) . 

REDUCED TRIDIAG^NAL MATRIX ELEMENTS Rt^W * 

CiFF DIAGtZiNAL - ** 

DIAGONAL »= *** 

* - The row number of the (reduced) tridiagonal matrix. 

♦♦ - Value of ^ne off-diagonal element for that row. 


*** “Value of the diagonal element for that row. 


Following the i>rinting of several linos containing reorthogonalization information and 
reduced tridiagonal matrix elements, when the FEER process has finished its computa- 
tions for the current point of interest, the header (see above) is printed once again, 
followed by a table which summari/.es all the eigensolutions found by FEER, This table 
, has seven columns, as follows; 

(1) Solution number. Tliis is simply a positive integer 1,2,3,... . 

(2) Order of extraction. These nmnbers indicate the order in which the 
tridiagonal matrix was constixicted. 

(3) Distance from center. This is the distance from the extracted eigenvalue 
to the neighborhood center (which is printed above the table) in the complex 
plane. The tabular values are sorted according to increasing distance from 
the center. 

(4) Real part of the extracted eigenvalue. 
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(5) Imaginary pail of tho oxtracled oIgcK^ ^^luo. 

(C) Thcorolicnl error oslimalo. This value must bo smnllor than tho Accuracy 
Criterion (sco above) for tlio olgcnsolution to bo accoptablo, 

(7) Status. A single word, *'acccpt" or "reject", to indicate tho result of tho 
accuracy test, A minus sign (-) is added to "reject" so that tho eye con 
• more rapidly distingulsli between ihe two words. 

Finally, this table is printed a second time, but with tiio rojected cigensolutlons deleted. 

For very small problems, tliero Is a Very Detailed Printout (VDP) option. This 
option was originally used to debug the complex FEEIl logic, and Is no longer required. 
DIAG 12 must be specified in tho executive control deck to Invoke the VDP option, and 
furthoi’moro, field G of a given EIGC bulk data continuation card must have a (floating 
point) value equal to or greater than the size of tlie stiffness matrix. Thus, field G of 
each EIGC continuation card should ordinarily be left blank, so that the VDP option will 
be suppressed. The actual printed output of this option consists of all vectors for each 
step of the complex PEER process, which is too extensive and detailed for normal usa’ 
purposes. 






OY 
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